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FOREWORD 


This  report  is  one  of  the  series  of  final  reports  on 
studies  carried  out  under  Contract  No.  SD-71  (ARP A). 

The  principal  contributor  to  this  paper  is  Richard  Kauth. 

This  report  was  illustrated  by  George  Zacharias  and 
prepared  for  press  by  Joyce  A.  Wegner. 

Respectfully  submitted, 
Laboratories  for  Applied  Sciences 
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1.  INTRODUCTION 


Many  of  the  proposed  or  existing  schemes  for  discrimina¬ 
tion  of  targets  against  background  depend  upon  knowledge  of  the  joint 
probability  distribution  of  one  or  more  parameters.  The  parameters  may 
have  the  units  of  time  (e.  g.  ,  crossing  time  distributions),  space 
(e.  g.  ,  position  as  in  the  neighborhood  modification  process),  spacial 
frequency  (Weiner  spectra),  or  radiance  (point  count  statistics).  What¬ 
ever  the  set  of  parameters  used,  whether  composite  or  "pure",  if 
there  are  n  parameters  then  each  observation  on  an  object  is  represented 
by  a  point  in  an  n  dimensional  space  of  the  parameters.  The  density 
of  such  points  as  a  function  of  the  parameters  is  the  joint  probability 
density  function  of  an  object  or  class  of  objects.  The  standard 
techniques  of  category  recognition  then  are  used  to  discriminate  among 
classes  of  objects.  Estimates  of  the  density  functions  for  different 
categories  of  objects  can  be  gotten  from  direct  observation  on  known 
objects,  or  from  assumptions  about  the  laws  governing  the  behavior 
of  the  objects  in  the  presence  of  certain  input  conditions,  combined 
with  estimates  of  the  probability  density  functions  of  the  input  conditions. 

In  this  paper  we  are  considering  principally  the  latter 

*/  k/ 

approach  applied  to  the  apparent  radiance  — '  ,/f  ,  of  the  earth's  atmo¬ 
sphere  in  several  wavelength  bands. 


—  Strictly  the  term  ’’radiance1'  refers  to  the  self  emission  of 
an  object,  while  "brightness”  or  "apparent  radiance”  refer  to  the  sum 
of  radiation  from  all  sources,  reflected,  transmitted  and  self  emitted. 
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2.  PURPOSE 


It  is  the  purpose  of  this  paper  to  prepare  the  groundwork 
for  the  computation  of  the  joint  probability  density  function  of  apparent 
radiances  in  several  spectral  regions,  seen  looking  down  at  the  atmo¬ 
sphere  from  above. 

The  physical  situation  is  as  follows: 

Radiation  from  the  sun  passes  through  the  atmosphere 
and  is  absorbed  and  scattered.  It  strikes  a  cloud  top  and  is  diffusely 
reflected  back  up  through  the  atmosphere  to  an  observer.  Towards 
longer  wavelengths  emission  may  also  become  significant  and  this 
contribution  may  be  handled  by  a  separate  term  involving  a  diffuse 
reflector  behind  an  emitter. 

The  question  of  probabilities  arises  because  of  almost 
random  changes  in  the  atmosphere.  Thus  the  parameters  of  an  atmo¬ 
spheric  model  are  only  know  probabilistically  as  are  the  radiance 
values. 

3.  THE  MATHEMATICAL  FRAMEWORK 


The  mathematical  framework  may  be  specified  as  follows: 
Given,  a  collection  of  variates, 

|  ,  j  e  n  ,  n  =  ^1,  .  .  . ,  n  | 

and  a  collection  of  dependent  variables  y.  =  y.  |  x j  j.  ,  i  €  n  and  given 
P  i  x.  I  ,  the  joint  probability  density  function  of  the  x  then  we  seek 


to} 

toJ 


,  the  joint  probability  function  of  all  the  y.,  and  further  we 


seek  ^  y^  \  ,  the  joint  probability  density  function  of  the  mth  subset 
i€  nr 

of  the  y^,  men. 
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If  condition  c)  is  not  met,  then  it  is  still  possible  to 
break  the  space  of  the  x.  into  regions  in  which  a  one  to  one  transfor¬ 
mation  does  exist  and  add  the  contributions  to  p  <  y.  >  from  each 

Y  l  if 

region. 

In  general  there  may  be  more  or  less  x/  s  naturally 

available  than  there  are  y.’s.  For  instance  we  may  be  interested  in 

only  two  of  the  y.,  say  y  and  y  ,  but  still  believe  that  four  of  the  x. 

i  1  ^  J 

must  be  allowed  to  vary.  In  this  case  we  can  with  no  loss  introduce 
y^  =x3»  y4  =x4*  The  Jacobian  them  simplifies  to 


9  (xi,x2) 


since  and  x^  will  be  explicit  functions  only  of  y^  and  y^  respectively 


In  the  other  case,  in  which  the  number  of  y  * s  is  greater 
than  the  number  of  x.’s,  a  functional  relationship  will  exist  among  the 


y^'  s  such  that  f  =  f  to  0.  Mathematically  this  implies  that  one  can 
have  perfect  discrimination,  i.  e.  ,  if  a  set  of  y  /  s  are  observed  not 
satisfying  f  |  y.|  =0,  then  one  is  not  observing  the  object  to  which 
the  model  applies.  Physically  it  means  that  one  must  choose  a  more 
sophisticated  model,  i.  e.  ,  allow  variation  of  more  parameters,  in 
order  to  get  a  more  realistic  answer. 

4.  THE  PHYSICAL  FRAMEWORK 


We  wish  now  to  show  that  the  physical  situation  can  be 
formulated  in  just  this  way,  i.  e.  ,  that  the  apparent  radiance  in  spectral 
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(4) 


region 


V 

ion  , 


can  be  written  as: 


M  =  A/i  |  xk  | 


in  which  the  are  variates  for  which  Px  j 


can  be  measured  or 
estimated  in  some  way,  and  that  the  conditions  a),  b),  and  d)  above 
are  satisfied.  (It  is  of  course  also  desirable  to  keep  violations  of 
condition  c)  to  a  minimum). 


The  apparent  radiance  seen  by  the  observer  in  some 
spectral  region  i  is  given  by: 


S.  R.  T. 
111 


+  N. 
1 


(5) 


in  which, 

S  is  the  solar  irradiance  at  the  top  of  the  earth’s  atmo- 
i 

sphere, 

R.  is  the  cloud  diffuse  reflectance, 
i 

T  is  the  transmission  from  the  top  of  the  atmosphere, 

l 

down  to  the  cloudfand  back  to  the  observer,  and 

N  is  the  radiance  of  the  atmosphere  above  the  cloud, 
i 


All  of  the  quantities  except  will  be  dependent  upon  the 
directions  of  sun  and/ or  observer.  These  directions  are  knowable  and 
thus  do  not  enter  this  problem  explicitly. 

The  assumption  that  one  may  simply  multiply  and  T^ 
is  justified  by  the  observation  that  R  varies  only  slowly  with  wavelength. 

The  problem  breaks  naturally  into  several  areas: 

1.  Reflection  of  the  clouds. 

2.  Transmission  and  radiance  of  the  atmosphere  above 
the  clouds. 

3.  Variates  of  the  model  atmosphere  used. 
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A.  Reflection  of  the  Clouds 

Here  there  are  several  approaches  available.  One  may 

let  the  R.  become  some  of  the  variates  x,  .  Then  it  is  necessary  to  es- 
1  k 

timate  from  data  the  joint  probability  distribution  of  the  R^  and  any  other 
parameter  which  also  may  be  a  variate  for  the  transmission  and/or  the 
radiance.  Such  a  variate  would  be  h,  the  cloud  height,  which  will  also 
effect  the  transmission  and  radiance  of  the  atmosphere  above  the  cloud. 
Thus  one  would  estimate  p  (  |  R^  j. ,  h);  and  from  this  p(h)  and  p  (  |  R^||h), 
the  marginal  probability  of  h,  and  the  conditional  probability  of  M 
given  h,  are  available.  A  second  approach  is  to  represent  R^  as  a 
function  of  other  variables  which  then  become  variates,  x^. 

For  example 

Rj  =  g(p)  f  (K.t)  (6) 

in  which  p  is  a  n  roughness  coefficient1'  which  describes  the  random 

variation  of  the  cloud  surface  orientation,  f  (K.t)  is  a  function  of  the 

x 

optical  thickness  of  the  cloud,  K.t.  K.  is  the  extinction  coefficient 

11 

for  the  spectral  region  i  and  t  is  the  cloud  thickness  in  cm.  of  water, 
f  (K.,t)  can  be  extracted  from  the  work  of  C.  Bartky  ~  —  f  and  H.  Brown  — ^ 
of  these  laboratories.  It  will  go  asympotically  to  a  maximum  as  t  in- 
c  reases. 

On  the  assumption  that  the  surface  orientation  of  the 
cloud  varies  quite  rapidly  within  the  field  of  view  of  the  observer,  g(p) 
will  be  a  smooth  function  of  p. 

It  appears  that  Bartky1  s  work  with  plane  parallel  clouds 
might  be  extended  to  develop  an  approximate  model  for  g(p). 
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B,  Transmission  and  Radiance 


The  transmission  in  each  spectral  region,  T^,  is  in 
general  a  function  of  the  cloud  height  and  of  the  temperature  and  the 
total  and  partial  pressures  of  the  absorbing  gases  all  along  the  path. 

The  distribution  of  temperature  and  pressures  as  a  function  of  height 
constitutes  a  model  atmosphere  and  it  is  the  coefficients  of  such  a 
model  which,  in  the  most  general  case,  must  be  assigned  probability 
density  functions  and  allowed  to  vary. 

For  the  radiance,  N  ,  similar  statements  hold,  with 
i 

one  possible  addition.  Because  the  self  absorption  of  emitted  radiation 
is  very  strong,  it  is  a  good  assumption  that  but  little  if  any  self 
emitted  radiation  will  reach  a  cloud  and  be  reflected  back  to  an  observer. 
The  observer  will  simply  see  a  contribution  from  the  self  emission  of 
the  atmosphere  above  the  cloud. 


Now  the  apparent  radiance  in  the  spectral  region  i  may 


be  rewritten  as: 


in  which 


A  =  S.  .(p)f(K.t)  T.({a.},  h) 

tN*H  • 

|  a.  |  =  the  collection  of  the  coefficients  or  parameters  of 
a  model  atmosphere,  and  h  =  the  cloud  top  height. 


T.  is  clearly  a  one  to  one  unique  transformation  with  re- 

57  1 

spect  to  h,—  but  not  necessarily  with  respect  to  all  of  the  a.,  and  one 

J 

may  be  forced  to  break  the  computation  into  pieces  as  mentioned  on 
page  4.  In  either  case,  x.  may  be  written  formally  as: 


^=yK(P(  t,  h. 


LAS-TR- 199-49 


7 


A  variety  of  assumptions  are  available  which  simplify 
the  problem  of  estimating  p^.  For  example,  assume  that  the  distribu¬ 
tion  of  the  a^  is  independent  of  p,  t  or  h,  but  that  the  distributions  of 
p  and  t  depend  separately  on  h.  Thus, 

Px  {*k}  =  Pa  {aj}  •  Ph(h)  •  pt(t  1  h)  •  pp(p  '  h)  •  (a 

If  in  addition  the  collection  a.  has  been  reduced  to  its  minimum 

J  6/ 

size  by  recourse  to  common  causes  —  of  pressure  and  temperature 
dependencies  then  the  a.'s  may  be  assumed  independent  and, 


p  i  a. I  =  n  p(a.) 

a(  J|  j  J 


Depending  on  the  data  or  theory  available,  and  p^ 
may  be  written  as  they  appear  in  equation  8,  or  as 


pp(p’h) 


Pt(t,h) 


respectively. 


Finally,  if  the  direct  obserable  is  p  (  |r^  |  ,  h)  then 

Eqs.  7  and  8  may  be  reformulated  to  take  account  of  that  fact.  In  all 
cases  the  system  of  equations  will  satisfy  the  physical  requirements 
implied  by  conditions  a),  b) ,  and  d),  the  conditions  of  uniqueness 


and  continuity  on 


non-zero  and  finiteness  on  J. 


f  ,  of  finiteness  and  continuity  on  _ * 

1  9  x: 


and  of 


C.  Other  Modifications 

The  constant  assumption  of  this  section  has  been  that 
the  transmission  (and  also  the  radiance)  model  used  would  require 
knowledge  of  the  pressures  and  temperature  all  along  the  path. 

The  only  difficulty  is  that  no  exact  and  usable  model  for 
7/ 

either  T.  or  N.  exist.  Gates—  has  done  a  rigorous  computation  for 


LAS-TR- 199-49 


8 


T.  for  water  vapor  and  CO^  for  a  homogeneous  path.  He  has  also  com¬ 
puted  the  transmission  for  long  slant  paths  using  the  average  pressures 

and  temperature  in  the  path  and  achieved  excellent  agreement  with  the 

*/  8/ 

CARDE  —■  data.  On  this  basis  E.  Walbridge  —  of  these  laboratories 

has  used  Gates  program  for  transmission  and  Bartky’s  results  for 
plane  homogeneous  clouds  to  compute  the  apparent  radiance  seen  by 
an  observer  in  the  1.  9p  and  2.  7}jl  bands,  for  a  variety  of  cloud  heights 
and  sun  and  observer  angles.  His  average  pressure  and  temperature 
through  the  path  above  the  cloud  is  computed  from  the  ARDC  atmosphere. 
As  pointed  cut  on  page  4  this  approach  using  two  dependent  variables 
but  only  one  independent  variate,  h,  gives  rise  to  a  functional  de¬ 
pendence  between  the  dependent  variables,  K  9  and/K  Additional 
quantities  might  be  used  as  variates  from  Walbridge1  svwork,  for  instance 
the  sun  and  observer  angles.  Then  9’/^Z  7^  Can  comPu^ec^ 

for  a  given  p^(h)  anc^  a  gi-ven  scan  geometry.  If  however  one  wishes  to 
consider  these  quantities  (sun  and  observer  angles)  to  be  known,  then 
for  a  more  realistic  estimate  of  py  some  additional  atmospheric  para¬ 
meters  will  have  to  be  allowed  as  variates. 

Attempting  to  extend  Gate’s  approach  by  admitting  the 

average  temperature  and  pressure  as  variate s, using  the  techniques 

implied  by  Eq.  3 ,  would  multiply  an  already  long  computer  program. 

A  useful  compromise  might  be  to  use  the  average  pressures 

and  the  average  temperature  as  variates  in  conjunction  with  one  of  the 

9-12/ 

many  less  rigorous  models  available - x  .  This  would  give  a  reason¬ 

able  number  of  variates  (from  4  to  6  variates  depending  on  whether  p 
and  t,  were  included)  which  should  give  reasonable  estimates  for  a 
similar  order  joint  probability  density  function.  The  main  advantage  of 


-  Referred  to  by  Gates,  Reference  7. 
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using  one  of  the  less  rigorous  models  is  that  a  much  greater  part  of 
the  work  might  be  done  analytically  before  numerical  computations 
are  undertaken. 


5.  CONCLUSION 

This  paper  has  attempted  to  demonstrate  that  the  com¬ 
putation  of  joint  density  functions  of  apparent  radiance  (as  opposed  to 
their  direct  measurement)  is  not  an  unreasonable  task.  It  should  be 
pointed  out  that  even  if  direct  measurements  are  made,  these  are 
always  limited  in  extent  by  costs  and  therefore  serve  as  tests  of  the 
atmospheric  models  used. 
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